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and  parameters  of  current  experimental  interest,  the  correction  to  the  rotation  required  for 
stability  is  less  than  one  percent.  Calculations  are  presented  for  a model  linear  trajectory. 
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I.  Introduction 

Currently,  there  is  a great  deal  of  interest  in  the  compression  of  magnetic  fields  and  plas- 
mas by  means  of  the  controlled  implosion  of  highly  electrically  conducting  hollow  metallic 
cylinders  (liners).  The  liners  are  launched  by  means  of  either  high  explosives,1,2  electromag- 
netic inductive  coupling,3— 6 or  hydraulic  compression7.  The  source  energy  which  launches 
the  liner  is  converted  into  kinetic  energy  of  the  liner  and  eventually  into  internal  energy  of  the 
enclosed  magnetic  field-plasma  mix.  After  the  liner  is  stopped  by  the  payload  pressure  (at  the 
inner  turning  point),  it  rebounds  and  the  energy  is  returned  to  the  liner,  which  moves  radially 
outward.  Thus,  the  liner  moves  inward  from  the  launch  point  to  a minimum  radius  and  max- 
imum pressure  condition  at  the  inner  turning  point,  and  then  back  out  again  to  the  launch  ra- 
dius. 

Different  types  of  liner  materials  have  been  considered  and  launched  upon  purely  mag- 
netic field  payloads  as  part  of  the  NRL  Linus  program.  These  include  solid  A1  and  Cu  liners 
and  liquid  metal  liners  of  a sodium-potassium  (NaK)  eutectic  mixture78.  The  current 
configurations  utilize  liquid  metal  liners. 

For  the  early  liners8,  a capacitor  bank  was  used  as  an  energy  source  and  the  liner  had 
two  free  surfaces,  the  inside  compression  surface  and  an  outer  free  surface.  In  an  elementary 
fashion  (detailed  analytical  and  numerical  analysis9  substantiates  this),  it  can  be  shown  that 
the  outer  free  surface  is  Rayleigh-Taylor  unstable  during  the  launch  phase  of  the  liner.  The 
inner  free  surface  is  stable  during  this  phase.  However,  in  the  vicinity  of  the  inner  turning 
point,  the  opposite  is  true:  the  inner  free  compression  surface  is  Rayleigh-Taylor  unstable, 
while  the  outer  free  surface  is  stable. 

The  inner  free  surface  can  be  stabilized  by  liner  rotation  about  the  cylinder  axis,  which 
introduces  a centripetal  acceleration  opposing  R.  With  sufficient  rotation,  the  inner  surface  can 
be  stabilized  throughout  the  entire  liner  trajectory.  It  would  then  be  unstable  only  because  of 
coupling  to  the  outer  surface  during  the  launch  phase.  This,  however,  can  be  prevented  by  go- 
ing to  a "captive"  liner7  where  the  outer  surface  remains  in  contact  with  a set  of  axial  or  radial 
pistons  driven  by  high  pressure  gas.  This  combination  of  rotation  and  hydraulic  piston  implo- 
sion leads  to  liners  which  are  theoretically  stable  against  Rayleigh-Taylor  modes. 
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While  rotational  stabilization  of  the  inner  free  surface  has  been  demonstrated 
theoretically  8 for  an  ideal  inviscid  liner  and  established  experimentally,9  the  precise  role  of  real 
liner  fluid  viscosity  has  remained  an  open  question.  In  this  paper  we  study  the  effects  of  liner 
viscosity  in  the  realistic  limit  of  large  Reynolds  number  [experimental  Reynolds  numbers  Re 
of  the  order  of  2 x 104  are  typical].  We  find  that  viscous  effects  are  concentrated  in  a thin 
boundary  layer  at  the  free  surface  whose  thickness  is  of  the  order  of  Re~1^,2.  The  bulk  of  the 
liner  behaves  essentially  inviscidly.  Within  the  viscous  boundary  layer,  the  rotational  speed  of 
the  liner  at  the  free  surface  is  reduced  from  that  which  would  occur  inviscidly.  The  boundary 
layer  thickness  and  reduction  in  angular  velocity  are  determined  as  quadratures  of  the  basic 
state  trajectory,  R(t).  Numerical  results  are  obtained  for  a model  prescribed  radial  trajectory 
and  presented  as  a function  of  the  compression  ratio  and  Reynolds  number. 


II.  Analytical  Treatment 

The  liquid  liner  is  assumed  to  undergo  a cylindrical  flow  corresponding  to  the  motion  of 
an  infinite  hollow  cylinder  which  is  imploding  while  rotating  aboufits  axis.  Rotation  is  neces- 
sary both  for  the  initial  formation  of  the  liner  and  to  stabilize  against  the  Rayleigh-Taylor  in- 
stabilities near  the  inner  turning  point.  In  this  study  we  will  concentrate  on  the  inner  free  sur- 
face and,  without  any  essential  change  in  the  physics,  will  considet  the  liner  to  be  infinitely 
thick. 

The  following  assumptions  are  made: 

1.  The  liner  is  launched  from  a state  of  rigid  rotation  with  angular  velocity,  fl0,  at 
time  t = 0. 

2.  The  liner  is  incompressible  (constant  density  p). 

3.  The  liner  kinematic  viscosity  v is  constant. 

4.  The  flow  is  laminar. 


5.  The  flow  is  one  dimensional 

6.  Surface  tension  is  neglected. 


Incompressibility  of  the  liner  in  this  case  implies 


-T-  ( ru ) =0,  (1) 

dr 

where  u(r,t)  is  the  radial  velocity  component.  The  radial  component  of  the  Navier-Stokes 
equation  governs  the  radial  motion. 


du  du 

— 1-  u — 

dt  dr 
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p dr  dr  r dr 


where  v (r,i)  is  the  zonal  velocity  and  p(r,t)  is  the  pressure.  The  angular  momentum  equation, 
which  is  formed  from  the  zonal  component  of  the  Navier-Stokes  equation,  is  given  as 

dr\  d , , 9 19,,  ,,, 

(3) 

Substitution  of  Eq.  (1)  into  Eq.  (2)  reveals  that  the  right-hand  side  of  the  radial  momentum 
equation  vanishes.  Viscous  effects  arise  only  through  coupling  with  the  zonal  momentum 
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equation  (3)  by  virtue  of  the  centripetal  acceleration  term,  -v2/r.  An  additional  viscous  effect 
arises  at  the  free  surface,  where  the  normal  stress  is  balanced.  This  will  be  discussed  later  in 
•he  paper  when  we  consider  the  boundary  conditions. 

The  coupled  Eqs.  (1-3)  should  be  solved  consistent  with  proper  initial  and  boundary  con- 
ditions to  provide  u(r.t),  \ (r,t)  and  p(r,t).  In  order  to  emphasize  the  main  physics  of  this 
problem  in  its  most  simple  form,  we  will  instead  consider  the  radial  trajectory  of  the  liner  to  be 
prescribed,  and  generate  consistent  solutions  to  the  angular  momentum  equation.  The  viscous 
evolution  of  the  angular  momentum  distribution  is  direct,  while  the  viscous  effects  upon  the 
radial  velocity  field  are  mainly  a result  of  the  coupling  to  the  zonal  equation  through  the  cen- 
tripetal acceleration  term,  — (cv)V. 

A.  Radial  Trajectory 

We  therefore  consider  a prescribed  radial  trajectory  of  the  free  surface 

R = RU/t0\ «),  (4) 

where  r0,  the  hydrodynamic  time,  represents  the  throw  time  (time  required  to  go  from  the  ini- 
tial radius  R0  e 1/2  to  the  minimum  radius  R0 ),  and  < 1/2  represents  the  ratio  of  minimum  to 
maximum  radius  (compression  ratio).  These  dimensions  and  parameters  are  indicated  in  Fig.  1. 

Integrating  the  continuity  equation  (1)  we  obtain 

ru  = R dR/dt  =A  (r),  ($) 

i.e.,  ru  is  a constant  throughout  the  liner  at  any  given  time  t.  Therefore,  the  specification  of 
the  free  surface  trajectory,  Eq.  (4),  completely  defines  the  entire  radial  velocity  field. 

The  general  approach  applied  in  this  study  is  valid  for  an  arbitrary  trajectory.  For  this 
study,  however,  we  take  a "parabolic"  form  for  the  prescribed  trajectory 


R ^ 

~ -1  + 

1 

1 

*0 

€ 

'o 

It  can  be  shown  that  this  type  of  trajectory  represents  a realistic  limit  of  a thin  freely  launched 
liner,  A /R  < < 1,  where  A is  the  liner  thickness.  The  numerical  results  to  be  presented  in 
this  study  will  be  associated  with  this  trajectory.  For  strict  consistency  with  our  previous  ap- 
proximation that  the  outer  radius  is  at  infinity,  it  is  necessary  that  viscous  effects  be  confined  to 
a layer  of  thickness  8 < < A. 

B.  Initial  and  Boundary  Conditions 

We  desire  a solution  to  Eq.  (3)  with  the  prescribed  trajectory  given  by  Eq.  (6).  Since  the 
liner  is  launched  from  a rigid  rotation  state  with  angular  velocity  ft0,  the  initial  condition  on 
A ( r,l ) - rv  is 

A (r,  0)  = r2n0,  r ^ R0e~i/2.  (7) 

A suitable  boundary  condition  must  also  be  applied  at  the  free  surface,  r - rt(r).  The 
free  surface  is  an  interface  between  a liquid  and  a magnetic  field  or  magnetic  field-plasma  mix- 
ture. In  either  case,  there  should  be  negligible  tangential  shear  stress  at  the  free  surface, 
and  the  normal  stress  in  the  liquid  liner  should  be  in  balance  with  the  plasma  and/or  magnetic 
pressure  p,  of  the  interior  volume  there.  Therefore  at  the  interface, 


0; 


(8) 
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"nt, 


P 


= Pi- 


The  vanishing  of  the  tangential  shear  stress  prescribed  by  Eq.  (8)  implies  that  at  r 
angular  momentum  satisfies 


A = Cr 2. 


(9) 

R (/)  the 
(10) 


Thus  the  zonal  velocity  field  is  locally  that  of  a rigid  body  rotation  with  arbitrary  angular  velo- 
city given  by  the  integration  constant  C,  which  is  in  general  a function  of  time. 


We  notice  that  the  boundary  condition,  Eq.  (8),  is  satisfied  by  the  initial  condition,  Eq. 
(7),  which  specifies  that  the  liner  is  launched  from  a state  of  rigid  rotation. 

Equations  (7)  and  (8)  provide  one  initial  condition  and  one  boundary  condition  for  the 
angular  momentum  equation,  Eq.  (3).  An  additional  boundary  condition  is  required  since  we 
have  a second  order  equation  in  r.  This  second  boundary  condition  arises  as  a far-field  condi- 
tion in  the  liner  for  r > R.  To  see  this  more  clearly,  we  observe  that  the  inviscid  solution 
(v  = 0)  of  Eq.  (3),  subject  to  the  initial  condition,  Eq.  (7),  satisfies  the  full  viscous  equation, 
but  not  the  boundary  condition,  Eq.  (8).  We  therefore  anticipate  that  far  from  the  free  surface, 
r > > R,  the  liner  will  behave  inviscidly.  A viscous  boundary  layer  will  exist  at  the  free  sur- 
face w'  recovers  the  boundary  condition,  Eq.  (8),  at  the  free  surface  and  matches  with  the 
inv  non  in  the  interior  of  the  liner.  This  matching  condition  provides  the  second  re- 

qi  lary  condition. 


.ire  2 illustrates  the  nature  of  this  viscous  boundary  layer.  The  liner  is  shown 
launched  at  time  t =0  with  a solid-body  rotation  corresponding  to  that  of  the  apparatus,  fl0. 
The  liner  is  also  depicted  at  a later  time,  r,  after  it  has  imploded  to  a radius  R(t).  The  inviscid 
zonal  velocity  variation  appropriate  at  that  time  is  indicated  by  the  solid  line.  This  variation 
corresponds  to  conservation  of  specific  angular  momentum  rv  for  each  of  the  fluid  lamellae 
composing  the  liner.  The  viscous  layer  or  boundary  layer  correction  is  shown  in  the  region  8 
by  means  of  the  dashed  line.  The  viscous  solution  satisfies  the  stressless  condition  at  the  sur- 
face and  matches  with  the  inviscid  distribution  in  the  interior.  Viscosity  is  seen  to  slow  down 
the  liner  interface  zonal  velocity.  This  is  important,  as  we  rely  on  the  strong  zonal  motion  of 
the  inside  surface  in  the  vicinity  of  the  turning  point  to  stablize  against  Rayleigh-Taylor  insta- 
bilites. 


C.  Lagrangian  Variables 

Because  of  the  motion  of  the  free  surface  and  the  inherent  competition  between  convec- 
tion with  geometric  convergence  and  diffusion  displayed  by  Eq.  (3),  it  is  more  convenient  to 
transform  to  Lagrangian  variables.  We  therefore  take 

i = r2  - R2  (r),  (11) 

which  is  proportional  to  the  displaced  area  relative  to  the  moving  interface  and  is  conserved  for 
fluid  shells  in  an  incompressible  liner.  This  transformation  has  been  used  previously  by  the 
present  authors  in  a related  problem  ", 


The  angular  momentum  equation  for  4(£,/)  becomes 
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U 
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(15) 

Equation  (15)  requires  that  the  solution  match  with  the  inviscid  solution  far  into  the  liner. 
D.  Nondimensionalization 

The  dependent  and  independent  variables  are  normalized  as  follows: 

A = A/ Rq  (1q,  t = t/ to, 

R =R/R0.  i =*//?o2- 

Two  dimensionless  parameters  characterize  the  problem: 


(17) 

(18) 


R2( 0)  ’ 

Re  — Rg  /4Wq  . 

Re,  a Reynolds  number,  represents  the  ratio  of  the  viscous  diffusion  time  Rq/Av  to  the  hydro- 
dynamic  time  associated  with  the  prescribed  radial  motion,  /0.  Since  this  is  a boundary  layer 
problem,  we  anticipate  that  the  viscous  zone  will  scale  as  Re-1/2.  As  previously  defined,  t is 
an  inverse  compression  ratio  representing  the  ratio  of  the  final  compression  area  to  its  initial 
value  for  the  prescribed  trajectory.  In  terms  of  these  non-dimensional  variables,  the  relevant 
equations  and  initial  and  boundary  conditions  become 
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at  large  where 
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(20) 
(21) 

(22) 

(23) 


A solution  to  Eqs.  (19-23)  valid  throughout  Re  - e parameter  space  represents  a formidable 
task.  For  large  values  of  Re,  the  situation  is  significantly  simplified.  This  limit  is  also  of  practi- 
cal interest,  as  the  Reynolds  numbers  associated  with  magnetic  compression  using  NaK  liners 
are  of  the  order  of  2 x 104. 
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E.  Asymptotic  Theory  at  Large  Reynolds  Number 

We  concentrate  on  the  large-Reynolds-number  portion  of  the  Re  — e parameter  space. 

In  particular,  we  seek  the  behavior  of  Eq.  (19-23)  in  the  following  asymptotic  limit: 

Re_1/2  < < « < 1, 

(24) 

and  will  apply  the  methods  of  matched  asymptotic  analysis  to  Eqs.  (19-23),  taking  advantage 
of  the  small  parameter  Re-1^2.  Henceforth,  we  omit  the  tilde  from  dependent  and  indepen- 
dent variables,  with  the  understanding  that  all  variables  are  dimensionless. 

a.  Outer  Solution 

fi. 

We  expand  the  angular  momentum  A in  a regular  perturbation  expansion  in  terms  of  the 
small  parameter  Re_1,/2: 

A({,r;Re.e)  ~ A<0>((.r,()  +Re_1/2/4(1)  (f,/;e)  + •••  . 

(25) 

To  lowest  order  in  Re-1^2,  Eqs.  (19-23)  become 

1 * 

© 

II 

© 

£Z>  <ro 

(26) 

with 

A (0)(£,0)  = £ + — , 

(27) 

dA(0>  A(  0) 

m — (o.')  = — . 

R2(t) 

(28) 

and 

A(0)  (f.r)  ~ - 

€ 

at  large  £.  The  solution  to  Eq.  (26)  which  satisfies  Eqs.  (27)  and  (29)  is 

(29) 

r 

1 1 

a ,0>  (£,/;«)  +-. 

€ 

(30) 

This  solution  does  not,  however,  satisfy  the  required  stressless  boundary  condition,  Eq.  (28). 

We  therefore  have  the  ingredients  for  a classical  boundary  layer.  It  is  necessary  to  rescale  the 
independent  variable  f to  recover  the  missing  boundary  condition  and  resolve  this  layer. 

b.  Inner  Solution 

The  necessary  stretching  of  £ is  given  as 

rj  = Re1/2f, 

(31) 

and  the  dependent  variable  is  written  as 

A = A (t) ,r;e  ). 

(32) 

In  terms  of  Eqs.  (31)  and  (32),  the  governing  equations,  Eqs.  (19-23),  become 

Y - I*2  </>  +nRe-1'2]  d2/j; 

d'  dy2 

A{  y.0)  = — + TiRe-1/2; 

(33) 

f 

■ 

(34) 

6 

k: 

r~— 



(0,/) 

-Re"1'2  A I0'0 ; 

(35) 

dr) 

R2(/) 

A (r),t) 

+ TjRe-1'2,!!  > > 1 . 

(36) 

The  initial  and  loundary  conditions  satisfied  by  Eqs.  (34)  and  (36)  are  evidently  the 
inner  representation  of  the  inviscid  solution  and  satisfy  the  differential  equation  Eq.  (33).  It  is 
therefore  convenient  to  define  a new  angular  momentum  variable  by  means  of 

£ 

A(r),r,e)  represents  the  viscous  correction  to  the  solution.  Equations  (33-36)  become 

>2: 


A (T),/;«  ) = A ( ) — 


— + rjRe-1'2 


(37) 


with 


M = [R2  (r)  +T)Re-*/2]-^- 
dt  dr,2 


A (r) ,0 ) =0, 


— (0./)  = Re-1'2 
dr; 


A (0,r)  + — 
e 


R2U) 


- 1 


and 


lim  A ( rj.l ) = 0. 


(38) 

(39) 

(40) 

(41) 


We  now  expand  A as  a regular  perturbation  expansion  in  powers  of  the  small  parameter 
Re-1'2: 


A (r),t,Re,e  ) — A <0)  (r),/;e)  + Re  ,/2/f(1)  (r),r;e)  + 
To  lowest  order  in  Re-’/2  we  find  that 

J<0)  =0, 

and  to  next  order  we  obtain 

= R2  (/) 


d2A(l) 


with 


dA 


dr) 


d'  dr) 2 

A(t)  (tj.0)  =0, 

(0,r)  = 1/e  R2(t)  - 1. 


and 


lim  A (l  * (tj,/)  = 0, 

T)—<» 


Defining  a new  timelike  variable,  the  age, 

H = f'/?2  (!')£//'. 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 


reduces  Eq.  (44)  to  a diffusion  equation  which  is  identical  to  that  governing  diffusion  in  a slab. 
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HA 


(l ) 


27U) 


i)  A 


, - , (49) 

cb/ 2 

It  should  be  borne  in  mind  that  tj  is  a dimensionless  area  and  w is  a stretched  time  variable 
which  is  weighted  over  the  trajectory.  That  is,  even  though  Eq.  (49)  has  the  form  of  a slab 
diffusion,  it  still  includes  effects  of  geometric  convergence  consistent  with  the  liner  motion 
asymptotically  correctly  for  the  large  Reynolds  number  considered  here. 


The  solution  to  Eq.  (49)  with  the  conditions  of  Eqs.  (45-47)  is  written  as 


<1 ) 


(t),H;«) -j-  f 


1 "r  rfr  p~V  2/4  («  — r ) 


[1/e/?2  (t  ) - 1]. 


(50) 


'/"•*()  («  - T ) 1/2 

Equation  (50),  the  viscous  correction  to  the  angular  momentum,  is  given  as  a quadrature  over 
the  basic  state  trajectory  R (t  ).  That  is,  A 11  * is  a function  not  only  of  the  compression  ratio  e, 
but  of  the  details  of  the  full  trajectory  in  terms  of  the  age  defined  in  Eq.  (48). 


In  summary,  the  angular  momentum  distribution  is  from  Eqs.  (37)  and  (50) 

A (t),/; Re,t ) 


+ r,  Re-'/2 

6 


- Re 


-1/2 


±s 


(he 


1 4 («  -7  ) 


(1  /tR2  - 1 ). 


(51) 


s/rr  1)  (ft  - r ) 1/2 

Of  particular  interest  is  the  angular  momentum  at  the  free  surface.  This  is  basic  to  stabilization 
of  the  Rayleigh-Taylor  instability  which  could  occur  near  the  inner  turning  point  / = i0.  Tak- 
ing t|  — 0 we  find 


A (0,/;Re,t ) 


1 


- Re 


—12  1 


f 


(h 

1 

(«  -r)"2 

e/?2  (r  ) 

- 1 


Therefore,  the  angular  velocity  of  the  free  surface  is  given  as 


(1 


I 


R 2 (ft ) 


± - Re"1  2 


J 


(It 


Jo  U>  - t) 


I 2 


I 


e /?  2 ( r ) 


— I 


(52) 


(53) 


= 1 - Re"1  2 * 


n , 


II 

f 

1 , 

7)  (n  - t)  1/2 

t R2  (t) 

^1  — Re-1  2/(t,f  ) 


(54) 


where  = 12 0/e /? 2 is  the  inviscid  surface  angular  velocity.  The  reduction  in  angular  velo- 
city of  the  free  surface  from  its  inviscid  value  is  seen  to  scale  as  Re-1  2. 

III.  RESULTS 

In  this  section  we  discuss  some  of  the  results  obtained  for  this  application.  Particular  em- 
phasis will  be  placed  on  the  viscous  pertubation  of  the  free-surface  angular  velocity,  since  this 
quantity  is  most  important  in  stable  and  well-controlled  basic  state  motion. 

The  quadratures  required  to  evaluate  the  scaled  specific  angular  momentum  profile  (Eq. 
(51)1  were  performed  numerically.  The  integration  variable  was  taken  to  be  t instead  of  »,  us- 
ing the  transformation  </«  = <liUlli/(li)  and  evaluating  tlli/tli  by  Eq.  (48).  Dividing  the  range 
of  integration  into  104  equal-time  intervals  and  using  the  trapezoid  rule  yielded  results  accurate 
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to  ~ 1%-  Landen’s  transformation11  was  used  to  evaluate  the  incomplete  elliptic  integral 
which  results  when  parabolic  trajectories  (Eq.  (23)]  are  employed  in  Eq.  (54). 

In  Fig.  3, /as  defined  in  Eq.  (54),  the  viscous  correction  to  the  free  surface  angular  velo- 
city, is  plotted  as  a function  of  time  for  various  values  of  t.  The  viscous  reduction  in  surface 
angular  velocity  is  seen  to  evolve  from  zero  at  t - 0 and  increase  with  time,  slowly  at  first,  and 
then  rapidly  near  the  inner  turning  point,  / - 1.  The  maximum  value,  of  order  1 in  all  cases,  is 
obtained  shortly  after  the  turning  point  for  this  model  trajectory.  As  the  liner  moves  outward 
towards  its  maximum  radius  at  t - 2.  the  correction  decreases.  The  correction  does  not  vanish 
entirely,  however,  even  at  t - 2 when  the  liner  has  returned  to  its  initial  radius.  Figure  3 
shows  that  there  is  signficant  « dependence  in  the  effects  produced  by  finite  viscosity.  In  par- 
ticular, we  see  that /tends  to  decrease  as  «— 0.  It  can  be  shown  in  this  limit  that  for  t*  1, 

/(/;«)  ~«1/2.  (55) 

while 

/(1;«  ) — « ,/4.  (56) 

provided  inequality  (24)  continues  to  hold. 

At  turnaround,  where  rotational  stabilization  is  most  critical, 

lim  e-,/4/0;«)  = 2.44.  (57) 

«— o 

Therefore  from  Eq.  (54), 

(n,„v  - n)/nmv  = 2.44Re-1/2e1/4  = 2.44(Re)^/2.  (58) 

where 

(Re)p#  =R0Rm3J^'0  (59) 

is  an  effective  Reynolds  number  related  to  both  the  minimum  and  maximum  radii  of  the  tra- 
jectory. 

But  it  is  also  clear  that 

lim  f(t ;«  ) = 0.  (60) 

«— i 

The  behavior  in  both  limits  is  shown  in  Fig.  4a,  where  we  plot  the  surface  angular  velocity 
correction  at  l = 1 as  a function  of  t.  Figure  4b  shows  the  same  curve,  scaled  by  ' 4 to  ex- 
hibit the  asymptotic  dependence  of  Eq.  (56). 

For  liquid  alkali  metals,  a typical  viscosity  is  «>=1  centipose.  Taking  t0  = 10“’  sec,' 
Rk  ■-  1 cm  and  « = 10~3  yields  (Re)^/2  = 10-3.  Thus  the  relative  change  in  surface  angular 
velocity  due  Ij  viscous  effects  is  always  less  than  one  percent  in  cases  of  practical  interest. 

IV.  SUMMARY 

Viscous  effects  for  imploding  liners  at  large  ratios  of  viscous  diffusion  time  to  hydro- 
dynamic  time  (high  Reynolds  number)  are  confined  to  thin  boundary  layers  near  the  stressless 
free  surface.  The  boundary  layer  thickness  scales  as  Re_l/2,  with  the  interior  region  of  the 
liner  behaving  inviscidly.  The  resulting  problem  has  been  solved  by  the  methods  of  matched 
asymptotic  analysis.  A boundary  layer  structure  has  been  obtained  which  satisfies  equations 
identical  to  those  associated  with  one-dimensional  slab  diffusion  while  properly  including  the 
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effects  of  geometrical  convergence.  The  results  obtained  indicate  that  the  free  surface  angular 
momentum  is  reduced  from  its  inviscid  value  by  an  amount  proportional  to  Re-1^2. 
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Fig.  1 — Geometry  of  the  model  used  to  study  viscous  effects 


at  the  free  surface  of  a rotating  imploding  liner. 


» 


Fig.  2 — Schematic  representation  of  the  viscous 
boundary  layer  at  the  free  surface  of  the  liner, 
showing  the  reduction  in  surface  zonal  velocity 
resulting  from  the  condition  of  zero  tangential 
shear  stress. 

1 1 
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Fig.  3 — Coefficient  of  the  viscous  correction 
correction  to  surface  zonal  velocity  vs  time  for 
various  compression  ratios. 
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Fig.  4(a)  — Plot  of  the  viacoui  correction  coefficient 
[Eq.  (54)]  at  turnaround  /(I*)  ▼«  e. 


